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Abstract. We study the interface between Regge behavior and DGLAP evolution in a non-perturbative 
model for the nucleon structure function based on a multipole pomeron exchange. This model provides the 
input for a subsequent DGLAP evolution that we calculate numerically. The soft input and its evolution 
give a good fit to the experimental data in the whole available range of x and Q^. 


1 Introduction 

There is a growing consensus on that nucleon structure 
functions (SF) at small virtualities are Regge behaved 
and that at large they follow QCD evolution j^. The 
two rerimes are incompatible because the DGLAP evo¬ 
lution Q, in general, changes the functional form of any 
Regge behavior (Regge models compatible with QCD evo¬ 
lution in a limited range of will be mentioned below). 
What remains completely unclear is the border between 
the two. 

Phenomenologically, the Regge pole approach to deep 
inelastic scattering implies that the structure functions are 
sums of powers in x, modulus logarithmic terms, each with 
a dependent residue factor. The rapid rise in of the 
structure functions, observed at HERA, initially was con¬ 
sidered as a sign of departure from the standard Regge be¬ 
havior. The reason was that the data, when fitted by a sin¬ 
gle “Regge-pomeron” term ^ x~^, where A is the pomeron 
intercept minus one, show that A = dirfi/x Q^) 

definitely rises with Q^. If so, Regge factorisation should 
be broken, since the “effective” pomeron trajectory ap¬ 
pears to be dependent. 

The above categoric conclusion was however based on 
a parametrisation of the SF by a single pomeron term, e.g. 
~ x~^. It is evident, that additional, decreasing (“sublead¬ 
ing”) terms, typical of any Regge pole model, will modify 
the properties of the leading one (pomeron). Moreover, 
by increasing the number of the free parameters and with 
complicated residue functions one may fit the data for 
arbitrarily large values of in a factorised form, with¬ 


out introducing any dependent trajectory, although 
it makes little sense to reach good at any rate with an 
indefinite number of free parameters introduced without 
any physical justification. 

The (in)dependence of the Regge behavior (trajectory) 
from the external masses (virtualities) comes from the fol¬ 
lowing arguments. The Regge asymptotic behavior 

( o \ Q;(t) 

-) , ( 1 ) 

So 7 

where ^(t) is the signature factor, (3{t,Q^) is the residue 
function, a{t) is the Regge trajectory and s,t are the 
Mandelstam variables, comes from a Sommerfeld-Watson 
transform over the partial-wave amplitude (see for exam¬ 
ple i) 


a{l,t,Q'^) 


I — a{t) 


( 2 ) 


By Regge factorisation, dependence can be intro¬ 
duced in the residue, but not in the trajectory. 

In this paper we are mainly interested in the small- 
X or high-energy behavior of the structure functions and 
photon-proton total cross sections, therefore we concen¬ 
trate on the pomeron (in terms of the Regge pole model) 
or singlet component of the SF, although, clearly, good 
fits and the evolution equation require the whole kine- 
matical region in x to be properly covered. Model for the 
pomeron, or singlet SF, and non leading contributions, or 
non-singlet SF, will be presented in Sec. ||. 
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Even though the role of the QCD evolution was stud¬ 
ied in a great number of papers (see e.g. jj] and refer¬ 
ences therein) the border and the interface between the 
non-perturbative Regge dynamics and perturbative QCD 
evolution is still determined merely by trial and error. De¬ 
pending on the scope of the study, one may prefer either 
approximate solutions of the DGLAP equation [H 

or computer calculations, for which purpose efficient nu¬ 
merical solutions and relevant computer codes are avail¬ 
able Q. 

Explicit analytical solutions are attractive for their 
simplicity and physical transparency. In paper , for ex¬ 
ample, the transition due to the evolution from a logarith¬ 
mic rise in x (dipole model) to a power-like was shown 
explicitly by means of an approximate solution of the 
DGLAP equation. Another example, in the Regge frame¬ 
work, is given by the CKMT model |P,[T^. 

Of particular interest are self-consistent solutions of 
the DGLAP evolution equations [||, i.e. those whose func¬ 
tional form does not change under the evolution. The 
product of logarithms (in x and in Q^) proved to be par¬ 
ticularly stable. In ref. 0 this stability for the calculated 
dependent coefficients, appearing in front of the fac¬ 
tors logarithmic in x, was checked experimentally. The 
stability of the combination of the logarithms 

A' 2 (x, Q^) = a-I-m^, ^ = In ^ In — (3) 

was noticed also by Buchmiiller and Haidt jn| from phe¬ 
nomenological analyses of the data. Ref. Qpapproaches 
this problem in a different spirit. The inclusion of a ”hard- 
pomeron” term, with Regge intercept about 1.4, allows 
for a fit to structure functions with a data point 

near 1. Moreover this new term is self-consistent, in the 
sense that a phenomenological parametrization of its Q^- 
dependence, at small turns out to be in agreement 
with perturbative evolution. However, the resulting dy¬ 
namics here is different from the one necessary to describe 
hadron-hadron scattering. 

Aiming at an unbiased extraction of the dependent 
factors appearing in the structure functions, the authors 
in iQ have fitted them for a number of different models 
of the pomeron, each appended by subleading terms as 
well as a large-x factor. The numerical values of the fitted 
coefficients for a large number of fixed values of (Figs. 
1, 2, 3 in El, see also Fig. 1 in ref. |^) can be used 
as “experimental values”, independent of any evolution 
scheme. 

The aim of this paper is to check the onset of QCD 
evolution for the model of small-x SF that was discussed 
earlier in ^,|^. This model contains a small (minimal in 
our opinion) number of free parameters, still it is feasi¬ 
ble and, supplemented by an evolution scheme, may have 
numerous applications. 

Our approach concerns the class of phenomenological 
parametrisation of structure functions describing the ex¬ 
perimental data in the non-perturbative region; using this 
parametrisation as input to the DGLAP equations one 
obtains a description of the experiment in the whole kine¬ 


matic rang e P, N[l- This genuine evolution, in contrast to 
global fits does not use data at large to con¬ 

strain the fit already at small Q^. In comparison to the 
DL Q model, where the pomeron has a hard nature, we 
choose a multipole pomeron with unit intercept. 

In the present paper the evolution is calculated numer¬ 
ically by means of the codes developed and published by 
Miyama and Kumano Q. 


2 Multipole pomeron model 

We analyze a unit intercept Multipole pomeron, rather 
than a supercritical one as in the CKMT case ||^. Here 
Multipole pomeron (MP) means that the pomeron is a 
multipole instead of just a simple pole. The advantages 
of this approach, both for hadron-hadron and for lepton- 
hadron scattering, were tested and discussed in numerous 
papers on that subject (see and 

references therein). The main point is that rising cross sec¬ 
tions (and rising SF) can be produced with a unit pomeron 
intercept. The number and the relative weight of the con¬ 
tributing multipoles is a very important question. In QCD 
it was studied in ref. . The conclusion of both the the¬ 
ory Po l and phenomenology is that a limited (moreover, 
small!) number of multipoles is sufficient in the present en¬ 
ergy (or x,Q^) range. In fact, a dipole (logarithmic rise), 
eventually with a minor tripole contribution (squared log¬ 
arithm) accounts for the large part of the data (higher 
order poles do not manifest). 

The dipole (and tripole) pomeron are typically “soft” 
objects. For the total cross sections they give fits close 
to those of a “soft” supercritical, i.e. with q;(0) « 1.08, 
pomeron. However, with an extra variable - t in hadron- 
hadron processes or in lepton-hadron reactions - the 
differences become essential. The result of the analysis 
below, including QCD evolution, makes this comparison 
more complete. 

Both non-perturbative inputs are written as sums of 
singlet and non-singlet terms: 

F2ix,Q^)=Fsix,Q^)+FNsix,Q^). (4) 

We use the same non-singlet term taking it from ref. 
as 

where, contrary to the original paper [ p^ , we disregard 
the difference between the light quarks. 

The singlet component of the SF, corresponding to a 
multipole (single-|-double-|-triple) pomeron is a sum of log¬ 
arithms: 
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Table 1. Parameters of the MP model 


MP model 

values of the parameters 

A 

0.8286 

B 

-0.1602.10“^ 

C 

0.8208.10“^ 

a 

0.7951 

b 

1.263 

c 

4.430 

a 

2.0 (fixed) 

p 

2.0 (fixed) 

7 

0.7711 

OiR 

0.415 (fixed) 

D 

1.434 

X^/dof 

1.08 



Fig. 1. Total 7 p cross-section versus W in the MP model. 


Fs(x, Q^) = [a (^)“ + B (^)^og (^) + 

_ ( 7 ) 

The same type of singlet component was used in Q. 

The real photon-proton total cross-section has the fol¬ 
lowing form: 


a^°p(W^) = dTT^asM ■ (A + B ■ log{W^)+ , . 
C ■ log^(fTi2) + D ■ ^ ’ 


where — rrip. 

The parameters of i)-® were determined from a com¬ 
bined fit of for IT > 4 GeV and structure function in 
the range of O.OdSGeT^ < < 5GeV^. 

We performed these fits to the set of the experimen¬ 
tal data using eqs. The values of the fitted 

parameters and of are quoted in Table 1. 

We have a rather small (eight) number of free param¬ 
eters. In the MP model we fix the reggeon intercept as in 
the CKMT model. Additionally a and (3 were fixed be¬ 
ing not sensitive to the ht. The resulting fits are shown in 
Figs. 1 and 2. 

As already mentioned, we chose the starting point for 
the QCD evolution in such a way as to have best fits with 
the smallest number of the free parameters. By trial we 
found that the value Qq = 6 GeV^ is preferable. In the 
numerical calculation, we assume that the only difference 
between the sea-quark and gluon distributions is a factor 
G(1 — x)~'^ in accordance with the dimensional counting 
rules. The value of the constant G is 3.8, the number of 
flavors is 4, A = 200 MeV. In contrast to Q, in our case 
the ratio of the gluon distribution to the singlet quark 
distribution is smaller by a factor of two. 

Figs. 2 and 3 show the extrapolation of the depen¬ 
dent “soft” Regge input as well as the results of our nu¬ 
merical calculations of the DGLAP equation by means of 
the NLO brute-force method, developed by Miyama and 
Kumano in ref. |||. thus extended to the high values of Q^, 
measured at |^f7 The input coincides with the structure 


function itself, supplemented by the aforesaid gluon distri¬ 
bution, since the corrections implied by the MS scheme 
are small and we neglected them. Although we are pri¬ 
marily interested in the small-a: behavior, large-a; data, 
including those from fixed-target experiments pl| were 
also taken into account since they can influence the small- 
X behavior through integration in the evolution equation. 


3 Conclusions 

We show that the world’s data on F 2 , down to — 0, 
are reasonably well described within the Regge model we 
have considered. The model succeed in reproducing the 
data for all x and a large interval, until « 6 GeV^, 
whereupon DGLAP equations must be used in order to 
describe correctly the evolution. What happens at the 
border between Regge regime and QCD evolution is clar¬ 
ified in the last icon of Fig. 2 and the first three frames of 
Fig. 3. While the starting point for the DGLAP evolution 
has been chosen at Qq = 6 GeV^, the extrapolation of 
the “soft” Regge behavior gives a satisfactory ht of the 
data until Q^ ^ 10 GeV^. In Fig.3, one can see that the 
ht, based on eqs. deteriorates when Q^ increases 

and, at Q^ = 20 GeT^, it is excluded by the data. There 
is however a region, a Q^ interval, where the two regimes, 
Regge and perturbative QCD, are compatible. The com¬ 
parison with other models, that consider only the small-a: 
region shows that the region of compatibil¬ 

ity shrinks when all the x-range is taken into account. To 
what extent this depends on the specihc model we have 
chosen, is an interesting question and will be studied else¬ 
where. 

To summarize, we hnd that the model described in eq. 
(^, for singlet SF supplemented by a non-singlet com¬ 
ponent as in eq. (||), together with the DGLAP evolu¬ 
tion, provides a simple and economic solution that could 
be useful for further practical applications, for example in 
nuclear physics. In addition, we can conclude that the soft 
input and its evolution give a good fit to the experimental 
data in whole available range of the variables x and Q^. 
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Fig. 3. The same as in Fig. 2. 








